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Abstract 

We discuss a new class of identities between correlation functions 
which arise from a local 2Z2 invariance of the partition function of 
the q— state Potts model on general graphs or lattices. Their common 
feature is to relate the thermal operators of the Potts model to some 
topological properties of the Fortuin-Kasteleyn clusters. In particular 
it turns out that any even correlation function can be expressed in 
terms of observables which probe the linking properties of these clus- 
ters. This generalises a class of analogous relations recently found in 
the Ising model. 



1 Introduction 



One of the simplest and more studied models of statistical mechanics is 
the q— state Potts model ||, which is the basic system symmetric under 
the permutation group of q elements. It can be defined on any connected 
graph G associating with each vertex i = 1, 2, . . . , v(G) the spin variable 
<7j = 1, 2, . . . , q . Its partition function at temperature T = 1//3 is taken to 
be 

Z G = J2e-W (1) 
M 

where TL = — J2tij) $0-^ with (ij) ranging over pairs of adjacent vertices in 
G. 

A long-standing problem is how to characterise geometrically in this 
kind of models the fluctuations near a critical point. Contrary to the naive 
expectation, the clusters made of adjacent sites with aligned spins do not 
play an important role in this respect. A different definition of cluster was 
proposed for the Ising model || and generalised to the q— state Potts model 
0]. These clusters are defined as adjacent sites with the same spin connected 
by bonds with probability p = 1 — . Within such a definition, these 
clusters behave correctly at the critical point, in the sense that their radius 
and the density of the percolating cluster scale with the correct critical 
exponents. 

The partition function (||) can be rewritten in terms of these clusters 
using the Fortuin Kasteleyn representation: || 

Z G (q,v)= £ ^ G V (G '\ (2) 
G'CG 

where v = = — 1; the summation is over all spanning subgraphs of 
G, namely the subgraphs made with all the vertices of G i.e. with v(G') = 
v(G); e(G') is the number of edges of G' , called active bonds, and k(G') 
the number of connected components or Fortuin-Kasteleyn (FK) clusters. 
This formulation of the partition function, sometimes called a di-chromatic 
polynomial, enables one to generalise q from positive integers to real and 
complex values. In particular q = corresponds to the tree percolation 
problem and q = 1 is the random percolation problem. 

In the present work we establish some new exact topological properties 
of these clusters. They arise from a class of identities between correlation 
functions which are valid for any graph G and come from the invariance of 
the partition function (|l|) with respect to a local 7Zi2 transformation. Their 
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origin suggests to call them Ward identities (see Sect. 2). Their common 
feature is to relate the thermal operators of the Potts model to some linking 
properties of the FK clusters and generalise a class of analogous relations 
recently found in the Ising model 0. 

It turns out that any even correlation function can be expressed in terms 
of a new kind of observables, which depend only on the closed circuits of 
these clusters (see Sect. 3). This suggests new, powerful methods to evaluate 
even correlators which can be easily implemented in numerical experiments: 
let G\ i G*2 ; • • • ; Gta , • • • be a sequence of spanning subgraphs in statistical 
equilibrium with the partition function (^), generated by any updating al- 
gorithm. Deleting all the edges which do not belong to any circuit of 
produces a dramatic simplification of the configuration with no information 
loss on thermal properties of the system. In the simplified configurations 
the evaluation of the new kind of observables is a much easier task. This 
leads to improved estimators of the thermal operators of this model. 



2 Ward Identities 



An elementary derivation of the simplest of these identities is the following. 
Let I be any edge of G. Denote by G\ the spanning subgraphs containing 
£ and by Gj the others. Clearly the set {G'} of all the spanning subgraphs 
may be written as the sum of two disjoint subsets {G'} = {G^} U {Gj}. 
Denote by 5e the 2&2 transformation on Zc{q,v) which deletes the edge £ 
from any subgraph of type G\ and adds it to any subgraph of type GJ. 
This is a symmetry of the partition function, namely, 

5 e [Z G (q,v)} = Z G (q,v) , (3) 

since J2ccg = J2g + cg + ^g~cg an< ^ ^ exchanges the two subsets {G^} «-> 
{Gj}. On the contrary the two partial sums are not invariant, but their way 
of transforming under 5i can be explicitly evaluated. We have, of course, 



e(Gj 



e(Gf) T 1 



(4) 



The action of on the number of clusters k(G^) depends on a topological 
property of t: if it belongs to a circuit of Gt , k is kept invariant 



k(G+)\ = k(G+) ; 

we call black bonds the edges with this property 
any closed circuit of Gf we have 



k(Gj 



k(Gj) + 1 



(5) 

If £ do not belongs to 
(6) 
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and call it grey bond or bridge. Using (Q), (||) and @ we obtain [] 



G+CG 



.1 



E [-vr G+ W + ^(l-vr G+ (£))]^)^ 



e(G+)„fc(G+) 



G+CG 



(7) 



where ttqi{1) denotes a projector on the black bonds defined as follows 

^ rg\ _ f 1 if ^ is a black bond ^ 
1 otherwise 

For the other partial sum we can write, more simply, 



E • 

GJCG 



e(G-) k(GJ) 



G+CG 



(9) 



From (||), (|7|) and (||) one obtains immediately 



(10) 



where (p(£)) = £ G+ P(Gj), with P(G') = v< G ^q k{ - G '^ /Z G {q, > 
Summing over all the edges I = 1, 2 . . . n = e(G) yields 



n 



■(6) 



fill 



where (e) = J2c e(G')P(G') is the mean number of edges of G', called active 
bonds and similarly (b) = J2c b(G')P(G'), where fe(G') is the number of 
black bonds in the subgraph G' . This is the first Ward identity. Since (e) is 
related to the internal energy E according to 



E 



d log Z G 



v+l 



d(5 v 
we can rewrite the above identity in the form 



(e) , 



E 



V + 1 -e(G) + iq - 1){V + '\b) 



v + q 



v(v + q) 



(12) 



(13) 



1 This way of reasoning is very similar to the one used to describe the basic properties 
of the Tutte polynomial of a graph Q , which is a simple variant of the partition function 
(0): see also [B. 
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which gives the internal energy E in terms of black bonds. This illustrates 
a common feature of this kind of Ward identities: all the thermal properties 
of the model are encoded in the black bond configurations, as we shall see 
shortly. The fact that the energy counts this special kind of bonds has been 



already observed some time ago [11 



As a little diversion, let us remark that if G were a tree T, then b(T) = 
b(G') = 0. This leads to E = — ^|e(T), which is a well known result of the 
q— state Potts model on any tree with sub-exponential growth [1C]; it also 
holds in the symmetric phase of any Bethe lattice . 



It is also interesting to study Eq. (11) in the pure percolation limit 



q — ► 1. Taking this limit together with the Taylor expansion of (e) 



< e > = ^fr + (« - x ) K e k ) - < e > WW + • • • * ( 14 ) 



yields the curious identity 

[(e)(k)-{ek)] a=1 , (15) 



(g ~ b )q=l 

— ^qrj — - l\w \"-/ - \ cf vj ? =i 



relating the mean number of bridges to the connected correlator between 
edges and clusters. 

Other Ward identities can be simply obtained by applying the above 2&2 
transformation to more than one edge. In particular, the two-edge trans- 
formation leads to the following expression of the specific heat C = ^ in 
terms of black bonds 



§ + £ + _ («-i)(«H-ir f 2W + 

(3 2 e(G) v(v + q) I w 

(b) 



+ 



(g - 1) [<6 2 ) - <6} 2 - -<&»]} (16) 



with (h) = J2g' h(G')P(G'), where h(G') denotes the number of the cutting 
pairs of the spanning subgraph G' , i.e. the pairs of black bonds such that 
their deletion increases the number of clusters. 



3 Cutting the bridges 

We come now to the other purpose of this paper, namely to show that not 
only the internal energy or the specific heat, but also any correlation function 
involving an even number of sites my be evaluated in terms of black bond 
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configurations B(G'), i.e. the graphs generated by deletion of all the bridges 
in the spanning subgraphs. To this end it is convenient to slightly generalise 
the model in the following way: we assign to any oriented edge (ij) of G 
a twist, namely an integer = 0, ±1, . . . , ±(q — 1) and define a modified 
Hamiltonian 

= -^Sai^j+nj , (17) 
(ij) 

where the sum <jj + is taken modulo q. An edge (ij) with ^ is said 
to be twisted. In order to get the Fortuin Kasteleyn representation of the 
twisted partition function Zq,t = e _/3Wr , we can rewrite it as 

Z G,r = EIE 1 + ™Wr y ) , (v = ^ - 1) , (18) 

M (ij) 

and expand the product as a sum of terms corresponding to the possible 
choices of summand from each factor. Each edge of the graph G can appear 
with a factor 1 or v5 aitaj+Tij ; in the latter case the edge is said to be an active 
bond. In a term with e(G') active bonds we have a factor ■y e ( G ") together 
with a product of delta functions which forces to zero every configuration G' 
in which all the spins connected by the active bonds do not satisfy the con- 
straint 5 = 1. If there are no twists in G' this constraint implies that all the 
spins of each cluster are aligned and the total number of spin configurations 
in a situation in which G' is formed by k(G') clusters is q k ( G \ If, on the 
contrary, there are twisted edges on G', we get a non- vanishing contribution 
only if there are no frustrations, namely if and only if the algebraic sum of 
the twists along any circuit of G' is zero modulo q. In order to select these 
subgraphs let us introduce a projector w T (G') defined as 

w (G 1 ) — \ ^ ^ ^' ^ as no f rustrat i° ns (ig) 
1 otherwise. 

Note that this projector picks out a subgraph G' regardless the location of 
its bridges, thus it can only be a function of the black bond configurations 
B(G>): 

m T (G') = w T (B(G')) . (20) 
In terms of this projector we have 

Zc,r^~ mT = E ^ e(G V (G,) ; (21) 

{a} G'CG 
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this leads to the fundamental identity 



where the expectation value is taken with respect to the standard, untwisted 
Hamiltonian. This defines a (complete) set of physical observables, one for 
each choice of r, which depend only on the black bond configurations and 
probe the linking properties of the FK clusters. 

Thus the problem reduces to show that any even correlation function 
can be expressed in terms of suitable combinations of tu T 's. This can be 



done, at least in principle, by writing the expectation value (22) in terms of 
spin variables 

{Wt ) = ^> ^ = (23) 

and then expanding the exponential. As examples, it is easily verified that a 



single twisted edge generates the first Ward identity (11) and the two-edge 
twist yields eq. (|l6|) . 



Note that eq. ( |23|) has the typical form of the expectation value of cor- 



relators of disorder observables [12]. Measuring disorder observables is not 
an easy task: one needs expectation values of operators which are expo- 
nentially suppressed on dominant configurations of the statistical ensemble. 
In particular, the numerical evaluation of operators involving a large num- 
ber of twisted links represents a very hard sampling problem. Actually Eq. 
( f20|) yields a simple, effective way out. Indeed w T (B(G')) is an improved 
estimator of the same observable. The evaluation of w T on the black bond 
configurations is much less noisy than in the spin configurations, because 
in the former case w T can only take the values and 1, while in the latter 
case the exponential of spin variables has a much larger variance, especially 
for large distance correlators. As a matter of fact, powerful algorithms to 
evaluate these new estimators have been constructed for the surface tension 



in 3D Ising model |T3 
7L2 gauge model |14|] - 
is now straightforward 



and for the Wilson loops Polyakov correlators in 3D 
Their extension to a general q— state Potts model 
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